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Abstract 

Using  the  finite  element  method,  lateral  currents  in  two  different  configurations  for  current  distribution  mapping  in  polymer  electrolyte 
fuel  cells  (PEFCs)  were  simulated  and  the  impact  on  the  accuracy  of  measurement  was  analysed.  The  measurement  techniques  were  of  a 
conventional  and  a  newly  developed  type,  based  on  a  segmented  bipolar  plate  (BPP),  respectively,  a  non-segmented  bipolar  plate  combined 
with  a  printed  circuit  board  (PCB).  In  both  cases,  neither  the  membrane  electrode  assemblies  (MEAs)  nor  the  gas  diffusion  layers  (GDLs)  were 
segmented  and  local  currents  were  detected  passively.  The  resistance  of  the  measurement  circuit  and  the  current  density  gradient  between 
neighbouring  segments  were  found  to  be  the  major  parameters  causing  current  spreading.  As  expected,  a  significantly  higher  uncertainty 
of  measurement  could  be  observed  and  experimentally  verified  for  the  non-segmented  bipolar  plate.  However,  the  accuracy  is  increasing 
with  increasing  homogeneity  of  the  current  density  distribution.  Achieving  a  uniform  utilisation  of  the  active  area  is  a  major  task  in  fuel 
cell  development  and  approaching  this  objective  also  improves  the  quality  of  measurement.  Consequently,  the  application  of  non-segmented 
bipolar  plates  as  highly  flexible  and  practical  measurement  technique  is  a  suitable  option  for  current  distribution  mapping  in  technically 
relevant  single  cells  and  fuel  cell  stacks. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

In  order  to  achieve  a  better  understanding  of  the  electro¬ 
chemical  processes  in  polymer  electrolyte  fuel  cells  (PEFCs) 
and  to  investigate  local  effects  of  operating  conditions  and 
hardware  configuration  on  cell  performance,  a  variety  of  cur¬ 
rent  distribution  mapping  techniques  were  developed  and  re¬ 
ported  in  the  literature.  Depending  on  the  type  of  fuel  cell 
to  be  analysed  (laboratory  scale  single  cells,  cells  of  tech¬ 
nical  relevance  and  fuel  cell  stacks)  and  on  the  phenomena 
to  be  studied,  the  individual  approaches  differ  with  regard  to 
methodology,  spatial  resolution  and  complexity. 

By  far,  the  majority  of  techniques  used  today  are  of  an  in¬ 
vasive  nature,  i.e.  current  distribution  measurement  is  associ¬ 
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ated  with  constructional  modifications  of  the  PEFC  hardware. 
Segmented  bipolar  plates  (BPPs)  in  particular  are  primarily 
applied  to  laboratory  scale  single  cells  and  consist  of  elec¬ 
trically  insulated  graphite  or  metal  segments.  Measurement 
of  local  currents  is  typically  carried  out  outside  the  fuel  cell 
either  passively  by  ohmic  resistors  [1,2],  Hall-effect  sensors 
[3,4]  or  current  transformers  [5],  or  actively  by  multichannel- 
potentiostats  [6,7].  Alternatively,  a  segmented  BPP  with  in¬ 
tegrated  Hall-effect  sensors  for  in  situ  current  detection  was 
documented  by  Wieser  et  al.  [8]. 

For  practical  applications,  BPPs  are  usually  manufactured 
in  a  thin  shape  for  volume,  weight  and  cost  reduction.  Conse¬ 
quently,  any  segmentation  is  difficult  to  achieve  if  mechanical 
strength  must  be  guaranteed  at  the  same  time.  Assuming  that 
the  ratio  between  through-plane  and  in-plane  conductivity  of 
thin  BPPs  is  high,  however,  lateral  currents  can  be  neglected 
and  segmentation  can  be  omitted  [9,10].  Following  this  as- 
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Nomenclature 

E 

vector  of  electric  field  strength  (V  m-1) 

i 

relative  segment  current  (%) 

I 

electric  current  (A) 

j 

vector  of  electric  current  density  (Am-2) 

R 

electric  resistance  (£2) 

u 

uncertainty  (%) 

Greek  letters 

<5 

boundary  layer  thickness  (m) 

<P 

scalar  electric  potential  (V) 

a 

electric  conductivity  (scalar/tensor)  (S  m-1) 

Subscripts 

0 

boundary  condition 

BPP 

bipolar  plate 

c 

contact 

eff 

effective  tensor  component 

GDL 

gas  diffusion  layer 

M 

measurement  circuit 

meas 

measurement 

mod 

modelling/simulation 

n 

segment  index  (1 1. .  .33) 

S 

shunt 

x,  y,z 

x-,y~,  z-direction 

sumption,  an  approach  with  “quasi-segmented”  BPPs  was  de¬ 
veloped  and  applied  to  single  cell  direct  methanol  fuel  cells  by 
Geiger  [9],  consisting  of  a  thin  non-segmented  BPP,  metallic 
current  collectors  on  the  backside  and  external  sensors.  The 
gap,  respectively,  the  uncovered  area  between  neighbouring 
current  collectors  additionally  suppresses  lateral  currents.  An 
alternative  approach  labelled  the  semi- segmented  plate  prin¬ 
ciple  makes  use  of  a  segmented  graphite  plate  attached  to 
a  thin,  non-segmented  BPP.  Local  currents  are  either  mea¬ 
sured  externally  as  with  the  “quasi- segmented”  cell  [10]  or 
the  ohmic  voltage  drop  across  the  individual  segments  is 
detected  in  situ  and  the  respective  currents  are  calculated 
[11]. 

A  further  invasive  method  for  current  distribution  map¬ 
ping  is  the  replacement  of  a  BPP  by  a  printed  circuit  board 
(PCB)  containing  gas  channels  and  current  collectors  on  the 
ribs  between  channels.  In  the  case  of  a  simple  PCB,  local  cur¬ 
rents  are  conducted  to  the  backside  by  plated  through-holes 
and  fed  to  a  sensing  device,  e.g.  a  multichannel-potentiostat 
[12,13].  If  the  PCB  is  manufactured  in  a  multilayer  design, 
current  sensors  such  as  shunt  resistors  [14]  or  a  combination 
of  coils  [15]  can  be  integrated  into  the  board  and  connected 
by  circuit  paths  in  the  intermediate  layers.  Obviously,  this 
configuration  is  particularly  capable  for  current  distribution 
measurements  in  fuel  cell  stacks  because  a  compact  design 
and  flexible  arrangement  of  the  wiring  is  possible. 

A  non-invasive  method  for  the  measurement  of  current 
density  distributions  in  fuel  cells  is  given  by  magnetotomog¬ 


raphy.  The  approach  is  based  on  the  fact  that  an  electric  cur¬ 
rent  flowing  in  a  conductor,  e.g.  a  piece  of  wire,  generates  a 
circular  magnetic  field.  Applied  to  fuel  cells,  sensors  detect 
the  magnetic  field  around  a  single  cell  or  a  stack  and  the  un¬ 
derlying  current  density  distribution  is  calculated  with  the  aid 
of  numerical  techniques.  Thus,  current  density  measurement 
without  any  drawbacks  on  fuel  cell  operating  conditions  is 
possible.  However,  magnetotomography  is  as  yet  in  an  early 
development  stage  and  only  allows  for  relative  measurements 
in  steady  state  [16]. 

A  major  problem  associated  with  invasive  measurement 
techniques  is  the  occurrence  of  lateral  currents,  most  no¬ 
tably  if  non-segmented  components  are  used.  As  a  conse¬ 
quence,  a  blurring  of  the  current  density  sets  in,  causing  an 
increased  measurement  uncertainty.  In  other  words,  current 
spreading  in  non-segmented  BPPs  and  gas  diffusion  layers 
(GDLs)  causes  the  measured  current  density  distribution  to 
be  more  homogeneous  than  it  appears  at  the  electrode.  Within 
the  scope  of  this  paper,  two  different  approaches  for  current 
distribution  mapping  are  chosen  and  the  respective  cell  con¬ 
figurations  are  modelled.  With  the  aid  of  the  finite  element 
method,  lateral  currents  are  calculated  and  drawbacks  on  the 
accuracy  of  measurement  are  analysed. 

2.  Scope  of  finite  element  modelling 

For  most  of  the  invasive  current  distribution  mapping  tech¬ 
niques  mentioned  above,  it  is  assumed  that  the  ratio  between 
through-plane  and  in-plane  conductivity  of  the  electrode  and 
GDL  is  high  due  to  the  small  layer  thickness.  Thus,  the  flow 
direction  of  the  electric  current  generated  at  the  triple  phase 
boundary  will  primarily  be  perpendicular  to  the  plane  of  the 
membrane  electrode  assembly  (ME A)  and  a  segmentation 
of  the  electrode  and  GDL  as  proposed  in  several  studies,  e.g. 
[7,11,17],  is  not  necessary.  On  one  hand,  avoiding  segmen¬ 
tation  at  the  place  of  electrochemical  current  generation  is 
desirable  to  preserve  true  fuel  cell  operation  characteristics 
and  avoid  highly  individualized  specialty  MEAs  [6] .  On  the 
other  hand,  however,  for  a  correct  interpretation  of  measured 
current  density  distributions,  it  must  be  resolved  to  what 
extent  the  assumption  of  negligible  current  spreading  is 
valid,  respectively,  what  measurement  uncertainty  must  be 
tolerated. 

Based  on  experimental  and  model  results,  Natarajan  and 
Van  Nguyen  [17]  recommend  a  segmentation  of  both  an¬ 
odic  and  cathodic  current  collectors,  respectively,  BPPs  as 
well  as  the  MEA  and  GDL  to  suppress  lateral  currents  and 
to  eliminate  effects  of  contact  resistance  variations  on  cur¬ 
rent  density  distribution.  However,  the  fuel  cell  configuration 
used  by  the  authors  is  highly  specific  and  only  one  single 
channel  is  regarded.  In  the  case  of  a  real  two-dimensional 
measurement  with  high  resolution,  a  significant  proportion 
of  the  active  area  may  be  electrochemically  inert  through  the 
insertion  of  insulation  layers  into  the  MEA.  Alternatively, 
potentiostatic  operation  of  individual  segments  helps  to  sup- 
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press  lateral  currents  by  partially  avoiding  a  gradient  of  the 
electric  field  strength  between  neighbouring  segments.  For 
high  spatial  resolution  or  simultaneous  current  distribution 
measurements  at  several  single  cells  in  a  fuel  cell  stack, 
however,  the  implementation  effort  is  rather  high  and  ex¬ 
pensive  due  to  the  multichannel-potentiostat  required  for  this 
approach.  Thus,  the  focus  of  this  study  is  concentrated  on 
non-segmented  MEAs  and  GDLs  and  on  passive  current  de¬ 
tection  techniques. 

Despite  of  an  inherently  low  measurement  resistance 
associated  with  the  magnetic  loop  current  detection  system 
described  by  Wieser  et  al.  [8],  a  maximum  measurement 
error  of  5%  was  reported.  Noponen  et  al.  [18]  used  a 
segmented  bipolar  plate  and  100  m^  precision  resistors  for 
current  distribution  measurements,  resulting  in  a  maximum 
uncertainty  around  10%.  From  these  examples  with  seg¬ 
mented  BPPs  and  non-segmented  MEAs/GDLs  documented 
in  the  literature,  it  can  be  concluded  that  the  occurrence 
of  lateral  currents  is  strongly  dependent  on  the  geometric 
configuration,  material  parameters  and  the  resistance  of  the 
measurement  circuit.  Furthermore,  it  is  obvious  from  these 
observations  that  significant  lateral  currents  can  be  expected 
in  the  case  of  non-segmented  BPPs.  The  thickness  of 
practical  BPPs  varies  in  the  range  of  0.1  mm  for  corrugated 
stainless  steel  sheeting  up  to  several  millimetres  for  graphite 
composite  plates.  Considering  not  only  the  geometric  but 
also  the  electrical  properties,  especially  graphite  BPPs  tend 
to  have  higher  intrinsic  conductivities  in  in-plane  direction 
than  through-plane.  Therefore,  depending  on  the  ratio 
between  through-plane  and  in-plane  conductivity  and  on 
the  geometric  parameters,  more  or  less  pronounced  current 
spreading  will  set  in.  Even  though  non-segmented,  respec¬ 
tively,  semi-segmented  BPPs  are  already  used  for  current 
distribution  measurement  in  practice,  only  rough  estimates 
concerning  lateral  currents  and  the  resulting  measurement 
uncertainties  have  been  published  for  these  methodologies 
[9]. 

For  a  quantitative  analysis  of  lateral  currents  in  different 
hardware  configurations,  two  techniques  have  been  modelled 
in  this  study.  The  first  one  depicted  in  Fig.  la  is  representa¬ 
tive  for  state  of  the  art  technologies  comprising  a  segmented 
BPP  at  the  anode  side  of  a  single  cell  with  individual  current 
collectors  and  external  sensors.  The  current  sensors  are  mod¬ 
elled  as  ohmic  resistors  Rm  and  neither  the  electrodes  nor  the 
GDL  are  segmented.  The  second  configuration  depicted  in 
Fig.  lb  represents  a  new  development  on  the  basis  of  preced¬ 
ing  approaches  for  current  distribution  mapping,  respectively, 
the  PCB  technology  and  the  semi- segmented  plate  principle. 
The  measurement  system  comprises  a  non-segmented  BPP  in 
combination  with  a  PCB  at  the  anode  side  of  a  PEFC.  Again, 
electrodes  and  GDLs  are  not  segmented.  Local  currents  are 
measured  by  means  of  shunt  resistors  Rs  integrated  into  the 
PCB  using  a  multilayer  design.  With  regard  to  stack  integra¬ 
tion,  the  PCB  contains  channels  for  a  coolant.  The  individual 
current  collectors  are  placed  on  the  lands  between  cooling 
channels  and  are  plated  with  gold  in  order  to  avoid  corrosion 


Fig.  1.  Schematic  of  cell  configurations  for  finite  element  modelling:  (a) 
segmented  bipolar  plate  and  (b)  non-segmented  bipolar  plate  and  printed 
circuit  board  at  the  anode  side  of  a  single  cell  (MEA,  membrane  electrode 
assembly;  GDL,  gas  diffusion  layer;  BPP,  bipolar  plate;  PCB,  printed  circuit 
board;  Rm,  measurement  circuit  resistor  and  Rs,  shunt  resistor). 


and  minimize  contact  resistance.  Depending  on  the  desired 
spatial  resolution,  several  current  collectors  can  be  grouped 
together  and  assigned  to  one  single  shunt  resistor.  Electri¬ 
cal  connections  between  current  collectors,  sensor  resistors 
and  a  conductive  layer  on  the  backside  of  the  PCB  are  pro¬ 
vided  by  plated  through-holes.  The  main  advantages  of  this 
newly  developed  method  can  be  summarized  as  follows:  (i) 
the  PCB  can  be  manufactured  in  a  shape  that  guarantees  for 
compatibility  with  the  anodic  bipolar  plate  in  a  single  cell, 
respectively,  with  adjacent  bipolar  plates  in  a  fuel  cell  stack; 
(ii)  the  sense  wires  associated  with  the  individual  resistors 
can  be  integrated  into  the  PCB  and  connected  to  a  data  ac¬ 
quisition  unit  from  the  edge  of  the  PCB;  (iii)  thus,  stack  in¬ 
tegration  at  any  position  is  possible;  (iv)  high  flexibility  due 
to  modular  configuration,  i.e.  different  BPPs  can  be  tested 
with  a  single  PCB  as  long  as  certain  design  parameters  (e.g. 
internal  manifolding)  remain  unchanged  and  (v)  as  with  a 
semi- segmented  plate,  the  fluid  dynamic  properties  for  the 
reactant  gases  as  well  as  the  electrical  and  thermal  conduc¬ 
tivities  of  the  PEFC  are  conserved  to  the  greatest  possible 
extent. 

Basis  for  the  simulations  was  a  PEFC  layout  provided 
by  the  Centre  for  Solar  Energy  and  Hydrogen  Research  in 
Ulm/Germany  with  100  cm2  active  area.  Anode  side  BPPs 
are  1 .5  mm  thick  with  a  serpentine  type  flowfield  as  shown  in 
Fig.  2a.  The  channel  width  is  1.0  mm,  channel  depth  0.3  mm 
and  the  width  of  the  landings  between  channels  is  1.0  mm. 
In  Fig.  2b,  the  parallel  channel  cooling  flowfield  integrated 
into  the  PCB  is  depicted.  Cooling  channels  are  1.5  mm  wide 
and  0.6  mm  deep;  the  width  of  landings  is  3.0  mm.  It  can 
also  be  seen  from  Fig.  2b  that  a  spatial  resolution  of  3  x  3 
segments  was  assumed,  i.e.  seven,  respectively,  eight  current 
collectors  are  grouped  together  and  connected  to  a  single 


174 


R.  Eckl  et  al.  /  Journal  of  Power  Sources  154  (2006)  171-179 


Fig.  2.  (a)  Anode  side  gas  flowfield  and  (b)  cooling  flowfield  of  printed 
circuit  board  with  current  collectors  configured  in  a  3  x  3  segment  pattern. 


shunt  resistor.1  In  an  analogous  manner,  the  segmented  BPP 
was  assumed  to  be  divided  into  3x3  segments  by  insula¬ 
tion  layers  (e.g.  an  epoxy  resin  material)  of  0.5  mm  thick¬ 
ness. 


1  It  should  be  noted  that  the  term  segment  is  used  in  this  study  both  for  seg¬ 
mented  as  well  as  for  non-segmented  BPPs.  In  the  case  of  a  non-segmented 
BPP,  the  area  beneath  one  or  more  current  collectors  associated  with  one 
single  current  sensor  is  considered  a  segment. 


3.  Model  equations,  methodology  and  material 
parameters 


The  current  generated  by  PEFC  is  DC  and  modelling  was 
carried  out  in  three  dimensions  for  stationary  states  in  conduc¬ 
tive  media.  Problems  of  this  type  are  described  by  Laplace’s 
equation 

V(trV^)  =  0,  (1) 


respectively, 


where  cp  represents  the  scalar  electric  potential  and  ax,  <jy, 

and  oz  are  the  components  of  the  electric  conductivity  tensor 

<t.  Furthermore,  Ohm’s  law  describes  the  correlation  between 

— ►  — ^ 

current  density  j  and  the  electric  field  strength  E : 

j  =  a  xE.  (3) 

The  correlation  between  electric  field  E  and  electric  potential 
cp  is: 


Thus,  follows  from  Eqs.  (3)  and  (4): 

j  =  -O  X  V<p,  (5) 


i.e.  the  simulation  of  electric  currents  requires  the  solution 
of  Laplace’s  equation.  For  this  purpose,  the  FEMLAB®  soft¬ 
ware  package  with  integrated  geometric  modelling  and  post¬ 
processing  capabilities  was  used  in  this  study. 

As  the  cathode  side  remains  unchanged  for  the  two  mea¬ 
surement  techniques,  the  simulation  can  be  restricted  to  the 
anode  side  alone.  The  objective  was  a  quantification  of  lat¬ 
eral  currents  in  the  components  GDL  and  BPP  and  no  elec¬ 
trochemical  processes  were  simulated.  Current  generation  at 
the  electrode  was  modelled  by  a  Neumann  type  boundary 
condition  with  a  defined  normal  component  of  the  inward 
current  density.  Furthermore,  a  uniform  compression  of  the 
GDL  was  assumed,  i.e.  no  distinction  was  made  between  ar¬ 
eas  beneath  channels  and  ribs  of  the  BPP.  From  this  follows  a 
constant  thickness  and  a  constant  through-plane  conductivity 
of  the  GDL  across  the  entire  active  area.  The  ohmic  resis¬ 
tances  of  the  individual  measurement  circuits  were  incorpo¬ 
rated  into  the  simulation  by  applying  Dirichlet  type  boundary 
conditions  at  the  current  collector  locations,  i.e.  by  defining 
equipotential  surfaces.  The  individual  voltage  values  were 
calculated  by  integrating  the  current  density  beneath  the  re¬ 
spective  current  collectors  and  multiplication  of  the  resulting 
segment  current  by  the  ohmic  resistance  of  the  measurement 
circuit.  As  the  numerical  solution  of  differential  equations  is 
carried  out  iteratively  in  FEMLAB®,  currents  and  potentials 
are  recalculated  within  each  iteration  step  until  convergence 
is  reached. 

The  contact  resistance  between  GDL  and  BPP  was 
modelled  by  thin  film  approximation  [19].  Assuming  a 
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negligible  tangential  component  of  the  current  density 
vector  and  a  constant  conductivity  ac  in  the  boundary 
layer,  a  one-dimensional  approximation  in  normal  direction 
(z-direction)  can  be  derived  from  Eq.  (2): 


-<*c 


This  equation  can  be  solved  analytically  by  applying  the 
potentials  <^bpp  and  <^gdl  to  the  upper  and  lower  edges  of 
the  boundary  layer  with  thickness  8: 


<^bpp  —  <^gdl  , 

- ; - z  +  <Pgdl- 


Together  with  Eq.  (5),  this  yields: 


<?BPP  -  <PG DL 

Jz  =  -°C - T - •  (8) 

For  the  simulation,  the  GDL  and  BPP  were  implemented 
in  two  separate  submodels  with  variables  <^gdl  and  <^bpp- 
Eq.  (8)  connects  the  submodels  by  defining  a  Neumann 
type  boundary  condition  in  each  case.  As  current  outflow 
from  the  GDL  represents  current  inflow  into  the  BPP,  the 
boundary  conditions  have  opposite  sign  in  the  respective 
submodels.  With  regard  to  the  surface  between  BPP  and 
current  collectors,  gold  plated  metal  sheets  with  a  negligible 
contribution  to  contact  resistance  were  assumed  in  the  sim¬ 
ulation  and  Dirichlet  type  boundary  conditions  as  described 
above  were  directly  applied  to  the  BPP  surface. 

Due  to  the  fact  that  each  segment  is  surrounded  by  at 
least  three  neighbouring  segments,  a  detailed  analysis  of 
lateral  currents  requires  three-dimensional  modelling.  For 
the  complex  flowfield  structure  of  the  BPP,  this  results  in 
considerable  computing  time  and  quickly  imposes  memory 
problems  when  a  standard  PC  is  used.  Furthermore,  the 
individual  components  feature  strongly  unequal  volume 
and  surface  dimensions  that  have  to  be  addressed  by  the 
discretisation  scheme.  An  approach  to  solve  these  problems 
is  the  numerical  volume  averaging  method  (NVAM),  i.e. 
replacing  details  of  the  original  structure  by  their  averaged 
counterparts  as  described  by  Roos  et  al.  [20].  For  this  pur¬ 
pose,  “virtual  experiments”  on  a  repetitive  region  are  carried 


out  and  tensor  components  for  a  representative  volume 
element  are  “measured”  using  the  finite  element  method. 

In  Fig.  3,  repetitive  and  representative  volume  elements  for 
the  given  anode  side  bipolar  plate  are  depicted.  The  repeti¬ 
tive  structure  comprises  a  channel  unit  and  half  a  rib  at  each 
side.  This  element  was  modelled  and  the  electric  conduc¬ 
tivities  of  the  BPP  material  adopted.  As  in  practical  exper¬ 
iments,  differential  voltages  were  applied  consecutively  in 
x-,  y-  and  z-direction  by  defining  Dirichlet  type  boundary 
conditions  at  opposite  surfaces.  At  the  remaining  surfaces, 
homogeneous  Neumann  conditions  with  zero  normal  cur¬ 
rent  density  were  applied.  The  current  flowing  through  the 
element  was  evaluated  numerically  in  each  direction  and 
finally  the  effective  conductivities  were  calculated  by  tak¬ 
ing  into  account  the  geometric  properties  of  the  represen¬ 
tative  volume  element.  Thus,  the  complex  geometry  of  the 
BPP  with  integrated  gas  channels  was  transformed  into  a 
simple  structure  that  could  be  simulated  with  a  standard 
desktop  computer.  Through  the  volume  averaging  process, 
all  the  information  about  the  underlying  BPP  structure  was 
cast  into  the  anisotropy  of  the  effective  conductivity  param¬ 
eters. 

In  the  PEFC  generation  used  for  this  study,  BPPs  are  man¬ 
ufactured  from  the  graphite-PVDF-compound  SIGRACET® 
BMA  5  with  specific  electrical  conductivities  of  100  S  cm-1 
in-plane  and  20Scm-1  through-plane.  NVAM  yields  for 
the  effective  conductivities  of  the  BPP  crx,eff  =  90  S  cm- 1 , 
Crieff  —  83  S  cm-1  and  crz  eff  =  16  S  cm-1.  Compared  to  the 
conductivity  in  y-direction,  the  higher  conductivity  in  x- 
direction  is  due  to  the  alignment  of  the  ribs.  As  GDL  material, 
SIGRACET®  GDL  10  BB  carbon  paper  was  assumed  with 
specific  electrical  in-plane  resistances  of  0.7  £2  square-1  in 
machine  direction  and  1.2  ^  square-1  in  cross  machine  di¬ 
rection.  Through-plane,  both  the  resistance  as  well  as  the 
thickness  of  the  GDL  are  strongly  dependent  on  the  mechan¬ 
ical  pressure.  For  the  simulation,  a  surface  pressure  of  17  bar 
was  assumed,  which  is  equivalent  to  a  compression  of  the 
GDL  to  75%  of  its  original  thickness  and  corresponds  to  the 
conditions  in  a  fuel  cell  stack.  From  this  follows  a  thickness 
of  0.32  mm  and  a  specific  electrical  through-plane  resistance 
of  9.1  m£2  cm2. 


Fig.  3.  Repetitive  geometric  structure  (left)  and  representative  volume  element  with  anisotropic  material  behaviour  (right)  for  the  anode  side  bipolar  plate 
(NVAM,  numerical  volume  averaging  method). 
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Table  1 

Electrical  conductivities  for  finite  element  modelling  of  lateral  cur¬ 
rents  (bipolar  plate  material  SIGRACET®  BMA  5,  gas  diffusion  layer 
SIGRACET®  GDL  10  BB,  mechanical  pressure  17  bar) 

Bipolar  plate  (S  cm-1) 

In-plane 

100 

Through-plane 

20 

Numerical  volume  averaging  (S  cm-1) 

x-Direction,  effective 

90 

y-Direction,  effective 

83 

z-Direction,  effective 

16 

Gas  diffusion  layer  (17  bar)  (S  cm-1) 

In-plane,  machine  direction 

45 

In-plane,  cross  machine  direction 

26 

Through-plane 

3.5 

Boundary  layer,  contact  resistance  (17  bar)  (S  cm  1 )  200 


Concerning  the  contact  resistance  between  GDL  and  BPP, 
no  values  were  available  for  the  combination  of  SIGRACET® 
GDL  10  BB  and  SIGRACET®  BMA  5.  For  this  reason, 
measured  data  for  graphite  and  carbon-fiber  paper  (Toray, 
3.5  wt%  Teflon)  as  documented  by  Mathias  et  al.  [21]  were 
applied.  At  a  surface  pressure  of  17  bar,  a  contact  resistance  of 
5.0  m£2  cm2  is  given.  This  corresponds  to  a  specific  electrical 
conductivity  of  200  S  cm-2,  which  represents  the  numerical 
value  of  crj&  in  Eq.  (8).  In  Table  1,  the  electrical  material  pa¬ 
rameters  for  the  individual  layers  used  in  the  simulation  are 
summarized  with  consistent  units.  BPP  and  GDL  material 
data  were  obtained  from  the  manufacturer  SGL  Technolo¬ 
gies  (Meitingen/Germany). 


4.  Simulation  results  for  the  segmented  bipolar  plate 
and  external  sensors 

For  the  purpose  of  a  worst-case  estimation,  a  strongly 
inhomogeneous  current  distribution  with  constant  current 
density  beneath  the  central  segment  and  zero  current  density 
beneath  the  surrounding  segments  was  assumed.  In  Fig.  4, 
the  simulation  results  expressed  as  relative  segment  currents 
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Fig.  4.  Relative  segment  currents  in  against  measurement  circuit  resistance 
Rm  and  position  of  segments  across  the  active  area  (segmented  bipolar  plate, 
3x3  segments). 


in  are  depicted  against  the  measurement  circuit  resistance 
Rm -  Relative  segment  currents  are  derived  from  the  integral 
segment  currents  In  through  the  equation  in  =  In/l o  x  100%, 
with  n  representing  the  segment  index  and  /o  representing 
the  total  current  flowing  into  the  structure,  i.e.  the  electro- 
chemically  generated  current  beneath  the  central  segment 
modelled  by  a  Neumann  type  boundary  condition  with 
constant  current  density  in  normal  direction. 

The  measurement  circuit  resistance  Rm  was  varied 
between  0  and  9  m£2,  a  range  that  reflects  many  state  of  the 
art  measurement  techniques.  It  can  be  noted  from  Fig.  4 
that  the  relative  segment  current  across  the  central  segment 
22  decreases  with  increasing  Rm  and  that  relative  currents 
across  the  surrounding  segments  increase.  Rm  =  0  can  be 
interpreted  as  galvanostatic  operation  mode  with  all  the 
individual  segment  voltages  forced  to  be  constant  by  a 
multichannel-potentiostat.  In  this  case,  about  7%  of  /o  are 
flowing  in  lateral  direction  within  the  GDL,  i.e.  the  measured 
current  density  distribution  is  always  more  homogeneous 
than  the  distribution  occurring  at  the  electrode  because  of 
the  inherent  through-plane  resistances  of  the  GDL  and  of 
the  graphite  segments.  Segments  12,  21,  23  and  32,  which 
are  conterminous  with  the  central  segment,  show  the  highest 
proportion  of  lateral  currents.  In  contrast,  lateral  currents 
across  the  corner  segments  11,  13,  31  and  33  are  almost  neg¬ 
ligible.  Furthermore,  it  can  be  noted  that  currents  occurring 
at  segments  21  and  23  are  higher  than  those  occurring  at 
segments  12  and  32.  This  is  due  to  the  fact  that  the  machine 
direction  of  the  GDL  with  higher  electric  conductivity  was 
applied  in  v-direction  of  the  model  geometry. 

The  strongly  inhomogeneous  current  density  distribution 
with  current  generation  only  beneath  the  central  segment  rep¬ 
resents  a  rather  theoretical  assumption.  Steep  gradients  can¬ 
not  be  excluded,  but  at  least  a  continuous  current  density 
characteristic  can  be  expected  in  reality.  Another  extreme 
example  is  defined  by  a  perfectly  homogeneous  distribution 
across  the  active  area.  Simulations  with  Rm  between  0  and 
9  mQ  revealed  no  significant  lateral  currents  for  this  operat¬ 
ing  state.  Thus,  it  can  be  concluded  that  current  spreading  is 
strongly  dependent  on  the  magnitude  of  Rm  and  on  the  cur¬ 
rent  gradient  between  neighbouring  segments,  respectively, 
the  gradient  within  a  domain  on  both  sides  of  the  insulation 
layer.  The  simulated  current  density  distributions — strongly 
inhomogeneous  and  perfectly  homogeneous — mark  the  up¬ 
per  and  lower  boundaries  for  the  measurement  uncertainty 
caused  by  lateral  currents  in  the  GDL  at  constant  Rm-  With 
increasing  current  density  gradient  and  consequently  with 
increasing  voltage  difference  between  neighbouring  current 
collector  surfaces,  the  measurement  error  is  also  increasing. 


5.  Simulation  results  for  the  non-segmented  bipolar 
plate  and  printed  circuit  board 

In  analogy  to  the  previous  chapter,  a  strongly  inhomo¬ 
geneous  current  density  distribution  with  current  generation 
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Fig.  5.  Relative  segment  currents  in  against  measurement  circuit  resistance 
Rm  and  position  of  segments  across  the  active  area  (non-segmented  bipolar 
plate  and  printed  circuit  board,  3x3  current  sensors). 

only  beneath  the  central  segment  was  modelled  for  a  first 
approach  and  the  simulation  results  are  expressed  as  rela¬ 
tive  segment  currents  in  Fig.  5.  It  should  be  noted  that  PCB 
measurement  circuit  resistances  include  the  individual  shunt 
resistors  Rs  for  current  sensing,  plated  through-holes  and 
internal  wiring.  Due  to  the  in-plane  conductivity  between 
neighbouring  segments  both  in  the  GDL  as  well  as  in  the 
non-segmented  BPP,  significantly  higher  lateral  currents  can 
be  observed.  The  highest  proportion  occurs  across  segments 
12, 21, 23  and  32  again.  However,  in  contrast  to  the  segmented 
BPP,  lateral  currents  across  the  corner  segments  11,  13,  31 
and  33  cannot  be  neglected  any  more  because  of  the  unhin¬ 
dered  cross-flow  of  electric  current  in  the  BPP.  Compared  to 
segments  12  and  32,  the  higher  currents  across  segments  21 
and  23  can  be  attributed  to  the  orientation  of  the  GDL  ma¬ 
chine  direction  and  BPP  channel  direction  along  the  v-axis 
of  the  model  geometry. 

In  order  to  test  the  methodology  with  non-segmented  BPP 
and  PCB  in  practice,  a  prototype  PCB  with  3x3  segments, 
cooling  flowfield  according  to  Fig.  2b  and  measurement  cir¬ 
cuit  resistances  of  7  mQ  (±0. 1  m£2)  was  manufactured.  With 
the  aid  of  this  prototype  PCB,  the  simulation  results  for  the 
non-segmented  PCB  can  be  verified  by  mounting  a  test  cell 
with  17  bar  clamping  pressure  and  defined  current  inflow  be¬ 
neath  the  central  segment  via  an  electrical  contact  sheet.  Re¬ 
sults  from  modelling  and  measurement  are  depicted  in  Fig.  6 
for  the  strongly  inhomogeneous  current  density  distribution, 
showing  a  good  accordance  between  simulation  and  exper¬ 
iment.  Currents  across  the  corner  segments  11,  13,  31  and 
33  are  slightly  overestimated  by  the  simulation  and  conse¬ 
quently  the  current  across  the  central  segment  is  calculated 
3.4%  lower  than  experimentally  observed.  However,  consid¬ 
ering  measurement  tolerances,  material  inhomogeneities  and 
a  scatter  of  material  parameters,  the  simulation  is  yielding 
correct  results  and  the  methodological  approach  developed 
in  this  study  is  applicable. 

As  with  the  segmented  BPP,  a  homogeneous  current  den¬ 
sity  distribution  was  modelled  for  the  combination  of  the 
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Fig.  6.  (a)  Simulation  and  (b)  measurement  of  relative  segment  currents  in 
%  for  strongly  inhomogeneous  current  density  distribution  (non-segmented 
bipolar  plate  and  printed  circuit  board,  3x3  current  sensors,  Rm  =  1  m£2). 

non-segmented  BPP  and  PCB.  The  results  from  the  previous 
chapter  can  be  assigned  without  any  constraints,  i.e.  no  lateral 
currents  can  be  observed.  Furthermore,  sensitivity  analysis 
revealed  that  the  current  density  gradient  between  neighbour¬ 
ing  segments  and  the  measurement  resistance  Rm  represent 
the  major  parameters  for  the  occurrence  of  current  spreading. 
For  practical  measurement  configurations  with  an  Rm  of  sev¬ 
eral  m£±  variations  of  GDL  thickness,  GDL  through-plane 
resistance  or  contact  resistances  in  a  simulated  ±20%  range 
have  a  minor  influence  because  of  the  strong  impact  of  Rm 
[22] .  This  behaviour  could  also  be  confirmed  experimentally 
by  varying  the  clamping  pressure  of  the  test  cell. 

6.  Comparison  of  measurement  techniques  and 
expansion  of  modelling  approach 

In  Fig.  7,  maximum  lateral  currents  for  the  two  ap¬ 
proaches,  i.e.  relative  currents  across  segments  21  and  23, 
are  depicted  against  Rm  for  the  strongly  inhomogeneous 
current  density  distribution.  Due  to  the  non-segmented  GDL, 
current  spreading  between  neighbouring  segments  can  be 
observed  for  both  measurement  techniques.  The  curve  for 
the  non-segmented  BPP  shows  a  relatively  steep  increase 
for  small  resistance  values,  followed  by  an  almost  constant 
characteristic.  A  significant  gain  in  measurement  accuracy 
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measurement  circuit  resistance  RV1  |mf2| 

Fig.  7.  Comparison  of  maximum  lateral  currents  between  neighbouring  seg¬ 
ments  for  the  worst-case  current  density  distribution  scenario  (BPP,  bipolar 
plate  and  PCB,  printed  circuit  board). 

is  only  possible  if  resistance  values  smaller  than  3  mQ  can 
be  realized  for  the  given  configuration.  The  curve  for  the 
segmented  BPP  in  contrast  is  monotonically  increasing 
within  the  considered  range  of  Rm  and  clearly  demonstrates 
the  inherently  higher  accuracy  of  this  approach. 

As  the  measurement  accuracy  strongly  depends  on  the 
current  density  gradient  between  neighbouring  segments 
and  current  distribution  mapping  only  yields  average 
values  for  individual  segments,  the  actual  uncertainty  of 
measurement  cannot  be  quantified  accurately.  Moreover,  the 
modelling  approach  developed  in  this  study  cannot  easily 
be  inverted,  i.e.  it  is  not  possible  to  evaluate  the  current 
density  distribution  at  the  electrode  in  a  one  step  procedure 
by  applying  the  measured  current  density  distribution  as 
Neumann  type  boundary  condition  at  the  current  collector 
surfaces.  For  the  electrode,  neither  Neumann  nor  Dirichlet 
type  boundary  conditions  can  be  defined  in  this  case  and 
thus  the  system  would  be  under-determined.  However,  it 
is  possible  to  estimate  uncertainties  for  a  detected  current 
density  distribution  by  applying  measured  current  density 
values  as  boundary  conditions  at  the  electrode  surface  and 
calculating  the  segment  currents  in  the  same  manner  as 
described  above  and  exemplarily  outlined  in  Fig.  8  [23]. 

Current  spreading  in  the  non-segmented  GDL,  and  as  the 
case  may  be  in  the  non-segmented  BPP,  results  in  a  ho¬ 
mogenisation  of  the  current  density  profile  that  can  be  quan¬ 
tified  by  a  measurement  uncertainty  un  =  (jmod,n  ~  7meas,rc) 
/ imesis,n  x  100%,  where  jmod,n  represents  the  simulated  aver¬ 
age  current  density  of  segment  n  at  the  current  collector  and 
7meas,n  represents  the  corresponding  average  current  density 
that  was  measured  and  defined  as  boundary  condition.  In  the 
example  according  to  Fig.  8,  an  experimental  result  obtained 
from  measurements  with  the  PCB  setup  is  analysed.  The  test 
cell  was  operated  in  galvanostatic  mode  with  H2  and  air  at  a 
total  cell  current  of  72.1  A.  The  ME  A  was  GORE  5620  with 
100  cm2  active  area,  the  other  PEFC  components  and  materi¬ 
als  were  the  same  as  described  in  the  Section  2  in  this  study. 
Maximum  relative  cross-flow  can  be  observed  towards  the 
bottom  left  segment  with  the  lowest  current  density,  originat¬ 
ing  from  the  surrounding  segments.  Local  currents  at  the  right 
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Fig.  8.  Approach  for  estimating  measurement  uncertainties  due  to  current 
spreading  in  non-segmented  gas  diffusion  layer  and  non-segmented  bipolar 
plate  for  a  measured  current  density  distribution  (non-segmented  bipolar 
plate  and  printed  circuit  board,  3x3  current  sensors). 


part  of  the  active  area  are  generally  measured  lower  than  they 
occur  in  reality  because  of  a  current  density  gradient  towards 
the  central  part  of  the  active  area.  Regarding  the  numerical 
values  for  the  measurement  uncertainties  in  Fig.  8,  it  must 
be  clearly  stated  that  these  can  only  give  estimates  because 
the  measured  current  density  distribution  serving  as  model 
input  already  represents  a  homogenised  version  of  the  true 
current  density  distribution.  Thus,  real  measurement  uncer¬ 
tainties  will  be  higher  than  those  specified  in  Fig.  8.  However, 
it  can  be  concluded  that  the  measurement  error  in  practical 
cases  is  significantly  lower  than  the  maximum  error  defined 
by  the  worst-case  current  density  distribution  scenario  ac¬ 
cording  to  Fig.  7.  Due  to  the  fact  that  increasing  the  unifor¬ 
mity  of  the  current  density  profile  obviously  coincides  with 
improving  the  quality  of  measurement,  the  application  of  a 
non-segmented  BPP  in  combination  with  a  PCB  can  be  re¬ 
garded  as  a  viable  diagnostic  tool. 


7.  Conclusions 

The  results  obtained  from  modelling  work  in  this  study 
clearly  demonstrate  the  impact  of  lateral  currents  on  the 
accuracy  of  measurement  for  different  current  distribution 
mapping  techniques.  Especially,  if  a  segmentation  of  the 
BPP  is  omitted,  significant  current  spreading  can  be  observed 
even  for  thin  BPPs.  With  the  aid  of  the  finite  element  method, 
a  detailed  three-dimensional  analysis  of  lateral  currents  is 
possible,  based  on  numerical  volume  averaging  of  complex 
flowfield  structures  and  on  thin  layer  approximation  of  con¬ 
tact  resistances.  It  could  be  demonstrated  experimentally  that 
the  modelling  approach  yields  correct  results  for  the  given 
geometries  and  material  parameters.  Assuming  a  strongly 
inhomogeneous  current  density  distribution  with  current 
generation  only  beneath  a  single  segment,  a  maximum 
measurement  uncertainty  can  be  defined  and  the  impact  of 
measurement  circuit  resistances  can  be  analysed.  For  the 
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techniques  chosen  in  this  study,  it  can  be  concluded  that 
the  accuracy  of  measurement  is  increasing  with  decreasing 
Rm  and  with  increasing  homogeneity  of  the  current  density 
distribution  generated  at  the  triple  phase  boundary. 

In  practice,  a  compromise  must  be  found  between  high  ac¬ 
curacy  and  technical  feasibility  of  a  measurement  technique. 
Segmenting  thin  BPPs  is  principally  possible,  however,  it  is  a 
challenging  and  tedious  task  to  guarantee  for  good  mechani¬ 
cal  stability  and  compatibility  to  an  existing  fuel  cell  layout, 
especially  if  higher  spatial  resolutions  than  3x3  segments 
shall  be  realized.  With  regard  to  stack  integration  and  flexi¬ 
bility,  the  approach  with  non-segmented  BPP  and  PCB  fea¬ 
tures  strong  benefits  compared  to  a  segmented  BPP.  Although 
the  measurement  uncertainty  associated  with  non-segmented 
BPPs  is  significantly  higher,  it  could  be  demonstrated  by  fi¬ 
nite  element  estimation  that  the  occurrence  of  lateral  currents 
is  less  pronounced  for  practical  current  density  distributions 
than  for  the  strongly  inhomogeneous  worst-case  scenario. 
Achieving  a  uniform  current  density  profile  under  a  variety 
of  operating  conditions  is  a  major  goal  in  developing  and 
constructing  PEFCs  and  coincides  with  a  gain  in  accuracy 
of  the  measurement  technique.  Thus,  the  combination  of  a 
non-segmented  BPP  and  PCB  in  fact  represents  a  viable  and 
flexible  alternative  to  existing  approaches  for  current  distribu¬ 
tion  mapping  in  technically  relevant  PEFCs  and  in  particular 
in  fuel  cell  stacks.  Irrespective  of  the  chosen  measurement 
technique,  an  estimation  of  current  spreading  is  possible  by 
means  of  finite  element  modelling  and  allows  for  a  correct 
interpretation  of  measured  current  density  distributions. 
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